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Pion-photon transition form factor using light-cone sum rules: 
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A global fit to the data from different collaborations (CELLO, CLEO, BaBar) on the pion- 
photon transition form factor is carried out using light-cone sum rules. The analysis includes the 
next-to-leading QCD radiative corrections and the twist-four contributions, while the main next-to- 
next-to-leading term and the twist-six contribution are taken into account in the form of theoretical 
uncertainties. We use the information extracted from the data to investigate the pivotal charac- 
teristics of the pion distribution amplitude. This is done by dividing the data into two sets: one 
containing all data up to 9 GeV^, whereas the other incorporates also the high-Q^ tail of the BaBar 
data. We find that it is not possible to accommodate into the fit these BaBar data points with the 
same accuracy and conclude that it is difficult to explain these data in the standard scheme of OCD. 
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FORM FACTOR F^ 



IN COLLINEAR QCD 



One of the most studied exelusive processes within 
QCD, based on collinear faetorization, is the pion-photon 
transition form factor with both photon virtualities being 
sufficiently large, see [l| for a review. The transition form 
factor is defined by the correlator of two electromagnetic 
currents 



d4ze-*«-^(^°(P) |T{j^(z)j,(0)}|0) 

= ie^,o.mi4F'^''''^{Q\<l^) (1) 



with = —q\ > 0, = — q| > 0, and can be reex- 
pressed in the form 

F^''^*"(Q2,g2)^^ / dxT{Q^,q'',fil,x) 



X v^^\x,t^l) + 0{6yQ') , (2) 



by virtue of collinear factorization, assuming that the 
photon momenta are sufficiently large Q^, 3> m^. Here 

N = \/2/3/x, /tt = 132 McV is the pion-dccay constant, 
and is the twist-four coupling. Then, the quark-gluon 
sub-processes, formulated in terms of the hard-scattering 
amplitude of twist-two, T, can be computed order-by- 
order of QCD perturbation theory: T = Tq + a^Ti + 
a^T2 + . . .. The radiative corrections in next-to- leading 
order (NLO), Ti, have been obtained in the /3o-part 
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of the contribution at the next-to-next-to-leading order 
level (NNLO/g), encoded in the amplitude T2, i.e., /3o-72/3, 
was calculated in 

The binding effects are separated out and absorbed 
into a universal pion distribution amplitude (DA) of 

twist-two, 'p^\x, fi^), defined Q by the matrix element^ 

(0|g(z)7.75[^,0]g(0)K(P)) ^ ^ iP^U Cdxe^^^^-P^ 

x^t'H^,A^F)- (3) 

The variation of ip^^\x,^p) with the factorization scale 
/ip is controlled by the Efremov-Radyushkin-Brodsky- 
Lepage (ERBL) evolution equation |2[; moreover the 
Gegenbauer harmonics {ipn{xY\ constitute the leading- 
order (LO) eigenfunctions of this equation. Therefore, it 
is useful to expand the pion DA in terms of these har- 
monics: 



^i2)(x,A.2)=^0(.T)+5^ 
' ^ r. 



an (iL^) 1pn{x) , (4) 



where ■>pn{x) =Gx{l~x)CfJ^ {2x~l) and ip^^{x)=')po{x) = 
6x(l — x) is the asymptotic pion DA 0- The nonpertur- 
bative information is contained in the coefficients a„ (/i^) 
with {n > 2) that have to be modeled or extracted form 
the data, including evolution effects to account for their 
/i^-dependcncc. They are usually reconstructed from 
the moments (^^)7r = Jo dx{2x — 1)^ ip'^\x, n"^) with 
(^°)7r = 1 that can be determined by employing, e.g., 
QCD sum rules (SR)s Q. We use here the pion DA 
proposed before in the framework of improved QCD SRs 
with nonlocal condensates (NLC-SRs) Q that yield a 



^ Gauge invariance is ensured by the longitudinal gauge link 
[z,0] = Pexp (ig Afi(T)dT'^) along a path-ordered lightlike 
contour. 
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"bunch" of admissible pion DAs with two harmonics that 
fix the coefficients 02 and 04. 



II. LIGHT-CONE SUM RULES FOR THE 
PROCESS 7*(Q^)7(q^ ~ 0) ^ 7r° 

The pion-photon transition involving two highly off- 
shell photons is not easily accessible to experiment. Ex- 
perimental information is mostly available for an asym- 
metric photon kinematics, with one of the photons having 
a virtuality close to zero — [sl-fioj. The calculation 
of this transition form factor within perturbative QCD is 
a precarious step because the quasi-real photon is emitted 
at large distances and has, therefore, a hadronic content 
calling for the application of nonperturbative techniques. 
An appropriate method is provided by light-cone sum 
rules (LCSRs) that supplements QCD perturbation 
theory with a dispersion relation for F'^ in the vari- 
able , taking then — ?> 0, whereas the large variable 
is kept fixed. Thus, one has 



s + q^ 
2 



(5) 



with the physical spectral density p{Q ,s) approaching 
at large s the perturbative one: 

1. 



-Im 



^7*7*^ (q2^-.s- ie) /N 



(6) 



Using quark-hadron duality, we obtain the following 
LCSR fnl: 



^2^7*7- (Q^)^9l /'^exp 



xN-p{Q\x)+ ^N-p{Q\x) (7) 
JO X 

with the spectral density p{Q^,x) = {Q^ + s)p^^{Q^,s), 
where s = xQ^/x and xq = Q^/{Q^ + so)- Note that the 
first term in ([7]) is associated with the hadronic content 
of a quasi- real photon at low s < sq, whereas the sec- 
ond term reproduces its point-like behavior at the higher 
value s > Sq. We adjust the hadronic threshold in the 
vector-meson channel to the value sq — 1.5 GcV^, us- 
ing nip = 0.77 GeV [l^l- We avoid to vary the Borel 
parameter in ([7)) and specify its value by virtue of 
= M|_p^/(a;)Q2 entering the two-point QCD sum 
rule for the p-meson with A/|_pt <E [0.5-^0.8] GeV^, where 
{x)q2 denotes some average value of x (at fixed Q^) in the 
integration region for the first integral on the right-hand 
side of Eq. © [ni[ll, i.e., xo{Q^) < {x}q2 < 1. 



III. MAIN INGREDIENTS OF THE LCSRS AND 
CONDITIONS OF THE DATA ANALYSIS 

It is convenient to invent for each term of the har- 
monics tpm a partial spectral density according to the 



definition ([6|) for the twist-two part [13], Pn\Q'^,s) = 
^ [(Tj «) ipn) {Q^, -s - ie)] . The general solution for 

Pn"^ in NLO was obtained in [l^ and corrected later 
(third line) in [ll]: 



+2v(n) In , „ 
\x Pf 



^ G„ii'i{x) +v{n)\ ^ &„iV'i(a 
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(=0,2, 



I '=0,1, 



,(8) 



with v{n), v^{n) being the eigenvalues of LO ERBL equa- 
tions, whereas Gni and bni are calculable triangular ma- 
trices (see [13, [3 for details). 

The inclusion of the NNLO^ contribution to the main 

partial spectral density /3o • p^^^ , derived from /3o • T2(5 [1] , 
was reahzed in It turns out that, taken together 

with the positive effect of a more realistic Breit-Wigner 
ansatz for the meson resonance [l^ instead of using a 
(5-function, i.e., 5{s — ■m'^)F'^ as in ([7]), it is negative 
and about -7% at small ^ 2 GeV^ , decreasing rapidly 
to -2.5% at > 6 GeV^. Here the results are expanded 
to include the first three harmonics. 

On the other hand, the twist-six contribution to F'' '''^ 
was recently computed in [l^ using ~ 1.5 GeV^ and 
found to be very small. Using instead the more moderate 
value 0.75 GeV^ [l3|, it turns out to have almost 

the same magnitude as the NNLO^ term, but with the 
opposite sign. 

As already mentioned, we use here the BMS bunch of 



pion DAs with the central point a^^^(lGeV^) 



0.20 



and 



,BMS 



(IGcV^) = -0.14 (termed the BMS model) 



These pion DAs have their endpoints at a; = and x = 1 
suppressed — even relative to the asymptotic pion DA — 
and are in good agreement with the CLEO data on the 
pion-photon transition form factor as well as with the 
data for other pion observables (l6l - [l8j . 

The key features of our data-analysis are the following: 
(i) The NLO radiative corrections in the spectral density 
are included via the corrected expression (|5]), emphasiz- 
ing that this error does not affect our previous results 
in [13, [l6l - [l8j . The so-called default renormalization- 
scale setting is adopted and, accordingly, the factoriza- 
tion and the renormalization scales have been identified 
with the large photon virtuality Q^. (ii) The twist-four 
contribution is taken into account using for the effec- 
tive twist-four DA the asyrnptotic form (pi'*-' (a;, p^) = 
(80/3)52(p2)a;2(i _ a.)2 ^ ^Iso admit a sig- 

nificant variation of the parameter = 0.19 GeV^ in 
the range 0.15 GeV^ to 0.23 GeV^, referring for details 
to [TtI , and taking into account its evolution with p,^. 

Using a nonasymptotic form for ip^^ would not change 
these results significantly [l^, [2^. (iii) The evolution 
effects of the coefficients a„ are also included in NLO, 
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employing the QCD scale parameters Aq^j-, = 370 MeV 
304 MeV, conforming with the NLO esti- 
0.118 [13. (iv) The NNLO/3 radiative 



and A'-qI^ = 

mate as(M|) ^ 

correction to the LCSR form factor [l^, d, Eq. 0, is 
incorporated together with the twist-six term, computed 
in [l5[, in terms of theoretical uncertainties. To be pre- 
cise, the calculation of the NNLO/3 term involves only the 
convolution of the hard-scattering amplitude Tp with the 
DA based on the three lowest harmonics. This treatment 
makes sense due to the fact that for the average value of 
A'P{Q'^) ~ 0.75 GeV^, these two contributions almost 
mutually cancel and the net result is small — see FiglTJ 
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FIG. 1: Twist-six contribution (upper solid line in red) and 
NNLO/3 contribution (lower solid line in green), obtained with 
the BMS model, and their sum (dashed blue line). 

In fact, it decreases with from -1-0.004 at = 
1 GeV^ — where the twist-six term dominates — down to 
-0.003 at Q2 ^ 40 GeV^— where the NNL0,3 cor- 
rection starts prevailing. This particular behavior ap- 
plies only to the moderate value of the Borel param- 
eter A'P = 0.75 GeV^ [l3l, while for the larger value 
= 1.5 ± 0.5 GeV^, used in the twist-six term 
would be much smaller and the net result would be ev- 
erywhere negative and almost constant: w —0.004. 



Gegenbauer harmonics necessary to model the pion DA. 
Second, we determine the fiducial regions of the corre- 
sponding coefficients a,i. Third, we relate these regions 
with the pion DA and its characteristics: profiles, deriva- 
tives at the origin, and its moments. Finally, we confront 
the obtained results with the data of set-1 and set-2. 



How many harmonics should be taken into 
account? 



To answer this question, we confront the dependence 
of the fit quality on the number of the parameters of the 
involved harmonics and the associated statistical errors. 
The statistical errors in the parameter determination in- 
crease with their number for statistical reasons, while the 
Xndi initially decreases. Therefore, in order to achieve an 
acceptable compromise, one should use the lowest ac- 
ceptable number of harmonics. The dependence of the 
goodness of fit, Xndi^ the number n of the involved 
harmonics for the two data sets, is presented in Fig. [5] 
The goodness of fit for set-1 is only slightly decreasing 
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FIG. 2; Dependence of the goodness of fit Xndf = X^/ndf 
(ndf — number of degrees of freedom) on the number n 
of Gegenbauer harmonics shown as histograms: set-1 — solid 
(blue) bars; set-2 — higher dashed (red) bars. 



IV. DATA ANALYSIS 

Here we overview our fit procedure of all available ex- 
perimental data on the pion-photon transition form fac- 
tor F'^ '^^ , within the framework of LCSRs, as worked 
out in [l3| . The main goal of the fit is to extract the pion 
DA — the main low-energy pion characteristic — best com- 
patible with all the data. This is done fitting the form 
factor by varying the pion DA in terms of the Gegenbauer 
coefficients a„. To reveal the particular role of the new 
high-Q^ BaBar data in the fit, we perform our analysis 
utilizing two different data sets. The first set (set-1) con- 
tains ah available data from CELLO Q, CLEG and 
BaBar [l3| that belong to the Q^-window [1 ^ 9] Gev^. 
The second set (set-2) comprises all data in the range 
[1 40] GeV^. First, we define the optimal number of 



with n and remains almost stable after n = 3. Thus, 2 
to 3 parameters are actually enough to describe all data 
in this region with Xndi ~ 0-5- In contrast, the data de- 
scription of set-2 is only possible with a Xndf value 2 or 3 
times larger — even if we include more harmonics. To fit 
all the data, we are forced to consider at least 3 param- 
eters with Xndf ~ 1- To have an even better description 
with a goodness of fit approximately equal to 0.8, we have 
to employ 4 parameters. Further increase of the number 
n will not provide any improvement. 

However, for the sake of comparison of the results, one 
should use the same fit model of pion DA, which means 
that the most appropriate number of harmonics may be 
fixed to 3. Best- fit curves for both data sets are shown in 
Fig. [3] as a bunch of form-factor predictions with errors 
stemming from the sum of the statistical error and the 
twist- four uncertainties. At high values of the momentum 
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transfer, the fit curve of the set-2 data — long dashed (red) 
hne — exceeds the 68% CL (confidential level) region of 
the set-1 data fit — solid (blue) line. This indicates that 
in the framework of LCSRs, the new BaBar data above 
9 GeV^ deviate from the low-Q^ data at the level of a la 
deviation and more. 
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FIG. 3: Best-fit curves to the experimental data for the tran- 
sition form factor in the framework of LCSRs: solid (blue) 
line — best-fit curve of set-1; strip bounded by dashed (blue) 
lines — 68% CL region; long dashed (red) line — best-fit curve 
of set-2; strip bounded by dashed dotted (red) lines — 68% CL 
region. Error bars show the sum of the statistical errors and 
the twist-four uncertainties. The experimental data are taken 
from the CELLO % (diamonds), CLEO (triangles), and 
BaBar [TO] (open crosses) experiments. 



B. Data analysis vs pion DA models 

Performing the data analysis, we obtain the best-fit 
values of the pion DA in the 68% CL region for a number 
of harmonics n = 2-;- 3. The 3D graphics of the confiden- 
tial regions for the 3 harmonics analysis were presented 
in our recent work in [l3| . whereas the best-fit values 
together with the statistical errors and the twist-four un- 
certainties are given in Table |I] below. We compare there 
our fit results with various pion DA models in terms of 
the goodness of fit Xndf analyzed sets of ex- 

perimental data. From the first two lines of Table HI we 
infer that the inclusion of the new high-Q^ BaBar data 
affect only the value of the parameter ae, while 02 and 
04 do not change significantly. Moreover, the good de- 
scription of the experimental data up to 9 GeV^ becomes 
appreciably worse after the inclusion of the high-Q^ tail 
of the BaBar data. The BMS pion DA stands out in the 
sense that it provides the best fit for set-1, while all other 
models cannot reproduce these data good enough. 

It is worth remarking that one of the models, obtained 
from fitting the experimental data within the Modified 
Factorization Scheme (MFS), has the value y^^f = 4.4 
in contrast to the result w 1.1 obtained in [24[. This 
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FIG. 4: (color online). Distorted \a error ellipses for set-1 
(upper panel) and set-2 (lower panel) from various exper- 
iments using different data-analysis procedures. 
These ellipses result from merging together the ellipses as- 
sociated with different values of the twist-four parameter in 
the range 5^ = 0.15-^0.23 GeV^. The slanted shaded (green) 
rectangle encloses the area of ai and 04 values determined 
by NLC-SRs [3, with the BMS pion DA being marked by K. 
The middle points of the ellipses (+ and T), the asymptotic 
DA (♦), the CZ DA (■), and Model III from [Q (A) are also 
marked. The range of values of 02, restricted by lattice sim- 
ulations, are indicated by vertical lines: [22| — dashed lines; 
[2^ 1 — dashed-dotted (blue) lines. All results are shown at the 
scale /^Iy = (2.4 GeV)'^, whereas the treatment of the Borel 
parameter M^(Q^) is explained in the text. Graphics taken 
from [H. 

discrepancy indicates that using the same pion DA in the 
framework of LCSRs and the MFS may lead for the same 
observable to incompatible results — a theoretical bias. 

Below, we consider in detail the results of the 2D anal- 
ysis in the (a2, 04) plane, presented in Fig. 21 with the up- 
per panel showing the results for set- 1, whereas the lower 
panel presents those for set-2. To this end, we calculate 
the Icr error ellipses^ by allowing the parameter to 
vary by 20% around the value 0.19 GeV^. The obtained 
error ellipses are then unified into a single (distorted) Ict 



^ We denote by a Icr ellipse (ellipsoid) a 68.27% confidence- level 
boundary. 
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TABLE I: Measures of goodness of fit of selected pion DA models (first column) with associated coefficients a„ (second column), 
used in the calculation of the pion-photon transition form factor by means of LCSRs. Note that the coefficients a„ are strongly 
correlated and the errors of a„ represent the maximal variation in the range of the fcr-region. The last two columns show the 
values of Xndf for the data in set-1 and for the whole set of the data, (set-2), respectively. All values of the coefficients a„ are 
given at the scale /isy ~ 2.4 GeV [25l ]. 



Model/Fit 


(a2, 04. . . .),,2_,,2 


Xndf, |1 - 9] GeV^ 


vLf, [1 - 40] GeV^ 

/vnaiJ L J 


3D fit, [i - 9] GeV^ 


(0.17 ± 0.11, -0.14 ± 0.18, 0.12 ± 0.17) 


0.4 




3D fit, [i - 40] GeV^ 


(0.18 ± 0.11, -0.17 ± 0.17, 0.31 ± 0.12) 




1.0 


NLG-SRs, BMS ]7} 


(0.141,-0.089) 


0.5 


3.1 


Model I [15] 


(0.084,0.137,0.088) 


> 2.8 


> 2.4 


Modif. fact. [24] 


(0.21,0.009) 


3.8 


4.4 


AdS/QCD, [26] 


(0.15,0.06,0.03,...) 


2.3 


2.8 


CZ [6] 


(0.394,0) 


32.3 


25.5 


Asympt. 


(0,0) 


4.7 


7.9 



ellipse shown in Fig. 21 To be specific, we consider the fol- 
lowing cases: (i) The result of combining the projections 
on the plane (02,04) of the 3D (3 parameter) data anal- 
ysis is represented by the largest ellipse — dashed (red) 
line with the middle point T. (ii) The analogous result 
of the 2D (2 parameter) data analysis in terms of 02 and 
04 is shown by the smaller ellipse (solid blue line) with 
the middle point + having the coordinates (0.15, —0.09) 
and Xndf ~ 0.5, that almost coincides with the middle 
point K of the parameter area determined by NLC-SRs 
(hi) The combination of the intersections with the 
(02, 04) plane of all 3D ellipsoids generated by the varia- 
tion around the central value of give rise to the smallest 
ellipse (thick line), entirely enclosed by the previous one. 

For convenience, the locations in the (02,04) plane of 
some characteristic pion DAs are also indicated in Fig. 
m These are the asymptotic DA (♦), the CZ model (■), 
and the projection of Model III from [l5j (a). Note that 
the slanted (green) rectangle, containing those values of 
02 and 04 that have been determined by NLC-SRs Q, is 
practically within both larger error ellipses and also over- 
lapping with the smallest one. Moreover, the BMS model 
DA >C stands out by lying inside of all la error ellipses. 
Thus, the theoretical predictions obtained from the 2D 
and the 3D data analyses conform with each other and 
agree at the level of Xndf — 0-5 with the results obtained 
from NLC-SRs Q- The calculated la error ellipses com- 
ply rather good with the boundaries for 02 extracted from 
two independent lattice simulations. The vertical dashed 
lines denote in both panels the older estimate from (2^ . 
while the very recent constraints from Ref. [l^] are rep- 
resented by the dashed-dotted (blue) vertical lines. 

From the lower panel of Fig. |4] it becomes evident that 
the situation changes significantly when including in the 
analysis the high-Q^ tail of the BaBar data In- 
deed, using the same designations as in the upper panel, 
we display the analogous unified error ellipses and ob- 
serve that the error ellipsoid has no intersection with 
the (02, 04) plane, whereas the composed error ellipse re- 
sulting from the 2D analysis (solid blue line) deviates 



from the region of negative values of 04 and moves inside 
its positive domain. At the same time, the fit quality 
deteriorates yielding Xndf ~ 2, as opposed to the value 
Xndf ~ 0.5 determined for set-1 of the data. As regards 
the unified la error ellipse of the 3D projections on the 
(02,04) plane (larger dashed red ellipse), its position re- 
mains unaffected, still enclosing most of the area of the 
02, 04 values computed with NLC-SRs — shaded (green) 
rectangle. 

The high quality of the data fit parallels the lattice 
findings, with the 3D error ellipse being almost entirely 
inside the boundaries from [22] (dashed vertical lines), 
while it also overlaps for the larger values of 02 with the 
range of values computed in [23] (dashed-dotted verti- 
cal lines). In contrast, the ellipse from the 2D analysis 
agrees very roughly with the small 02 window of [23| . 
sharing also only a small common area with the low end 
of the 02 region found in [2^ . Obviously, no agreement 
between the 2D and the 3D analysis is found. This dis- 
crepancy is also reflected in the values of the x-critcrion 
of the 2D fit model Xndf ^ 2 and that for the 3D model 
which turns out smaller by a factor of 2: Xndi ~ ^- This 
means that a pion DA, based only on 2 harmonics, is 
not sufficient to describe all the data on the pion-photon 
transition form factor. This deviating behavior of the re- 
sults, associated with the fits for set-1 and sct-2, shows 
up for a larger number of degrees of freedom, i.e., when 
including into the data analysis the next higher harmon- 
ics ipn with n = 6, 8, 10. But, for the case of the set-1 fit, 
this expansion does not improve any further the value 
Xndf — 0.5 — this remains approximately stable. More- 
over, the ItT-admissible regions in 2D, 3D, or 4D param- 
eterizations appear to be each embedded inside the other. 
In contrast, fitting the set-2 data, these new degrees of 
freedom lead to a decrease of Xndf, while the correspond- 
ing Icr-regions in the 2D, 3D, or 4D space, either do not 
overlap at all or intersect only marginally. 

It becomes obvious from Fig. 2] that Model III (a) 
from [Tsj has a projection on the (02,04) plane that lies 
outside of all considered la error ellipses of the data. 
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FIG. 5: Left. Comparison of the BMS pion DA bunch (shaded strip in green color) and of the BMS model (black solid line 
inside this strip) with the 3D fit to the experimental data on the pion-photon transition form factor. The solid blue line denotes 
the best-fit pion DA sample obtained from the analysis of set-1, with the dashed lines indicating the sum of the statistical 
errors of the fit and the twist-four uncertainties. Right. Analogous results obtained with set-2. 



However, seleeting for the Borel parameter the value 
A/^ = 1.5 GeV^, as in [l^, the agreement of this model 
with the data improves to the level of xidi ^ 



C. Pion DA characteristics 

The confidential region of the coefficients {a„}, ob- 
tained above, can be linked to any other characteristic of 
pion DA. The profiles of the pion DA ipTr{x), extracted 
in the 3D fit procedure, are shown in Fig. [SJ left panel — 
set-1; right panel — set-2. The BMS bunch (shaded green 
strip) and the BMS DA model (black solid line) are also 
shown in both cases. The inclusion into the data fit of the 
high-Q^ BaBar tail, causes a modification of the shape 
of the pion DA — sec Fig. [S] — giving support to our previ- 
ous observation that the BMS bunch is within the error 
range of the set-1 fit (left panel), while the best fit to 
set-2 differs considerably (right panel). In addition, the 
pion DA becomes endpoint enhanced, as opposed to the 
endpoint-suppressed BMS pion DA. The endpoint behav- 
ior can be characterized by its slope at the origin given by 
the derivative ^p'^{0) or, more adequately, by the so-called 
"integral derivative" D^'^'>ipTr{x), introduced in [23|- The 
integral derivative is the average derivative </9^(a;) defined 
by 




with the important property lim D^'^^ip[x) = </3^(0). 

Using a 3D confidential bound on the Gegenbauer co- 
efficients, we get the values of the derivatives V3^(0) and 
D^'^^Lp{QA), supplied in TablelUfor both data sets. These 
characteristics are shown together with the theoretical er- 
rors, the first being statistical and the second stemming 
from the twist-four uncertainty. We observe that these 



TABLE IL Comparison of the pion DA characteristics for 
the data of set-1 and set-2. 



data set 


[1 - 9] GeV^ 


[1 - 40] GeV^ 


v'Uo) 


20.2 ± 19.8 ± 1.1 


48.5 ± 11.4 ±0.4 


i3'2VT(0.4) 


6.6 ± 1.1 ±0.4 


8.1 ±0.7 ±0.3 


BMS DA 


Agreement 


No 


n 


2, 3 


3, 4 


Xndf 


0.53, 0.44 


1.0, 0.77 



characteristics can clearly differentiate the pion DAs gen- 
erated from set-1 and set-2. 



D. Combining Lattice constraints with 
experimental data 

The possibility to extract information on the mo- 
ment (^^)7r of the pion DA by combining lattice con- 
straints with the experimental data was first pointed 
out in and the following range of values was ex- 
tracted from the la error ellipse of the CLEO data @ 
in conjunction with the lattice constraints for (^^)^ from 
[H: {i"^)^ e [0.095 ± 0.134] at /x^t = 4 GeV^ and for 
= 0.7 GeV^. This procedure was refined in [11] in 
the following way: first, we expanded the result of the 
2D analysis to the ((^^)7r, (C^)7r) moments. Then, we de- 
termined the intersection of the confidential region (the 
area enclosed by a solid blue line) in Fig. |6]for set-1 (cf. 
Fig. SJ using the constraints from [23| and [231 ]. The 
intersection of these constraints, evaluated at the typi- 
cal lattice scale /.t^aj = 4 GeV^, and the experimental 
data leads to the following moment results, respectively, 

(i) (^2)^ £ [0.23 ± 0.29] and (^4)^ e [0.102 ^ 0.122], 

(ii) (^^)^ e [0.26 0.29] and (V),. G [0.11 ± 0.122]. 
These common validity ranges were extracted using a 
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1 
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0.20 0.22 0.24 0.26 0.28 0.30 0.32 

FIG. 6: (color online). Predictions for the moments (^^)7r 
and {£,'^)n at the lattice scale /^Lat ~ 

4 GeV^ The solid (blue) 
ellipse corresponds to our choice of M^, whereas the dashed 
(red) one results when using = 1.5 GeV^ The vertical 
lines show the range of values computed on the lattice: dashed 
line— [2^: dashed-dotted (violet) line— [111 • 

Q^-dependent Borel parameter — like everywhere in our 
analysis here and in [TjI- On the other hand, the value 
~ 1.5 GeV^ [1^, yields only a small intersection of 
the validity region extracted from set-1 (shown in F ig, p 
by the dashed red line) with the lattice constraints of [22] ■ 
This restricts the common region of validity to the value 
{£,^)Tr — 0.1, whereas there is no intersection at all with 
the lattice estimates from [2^. This obvious sensitivity 
of on the choice of the particular value of the Borel 
parameter gives additional support to our choice of 
the value of the Borel parameter. 

V. CONCLUSIONS 

We have presented here a global fit to the data on the 
pion-photon transition form factor, discussing further our 
recent analysis in . To get a precise measure of the 
influence of the high BaBar data on the form fac- 
tor and the pion DA, we divided the experimental data 
in two different sets with respect to Q^. Set 1 contains 
all data in the range [1 ~ 9] GeV^, whereas the second 
set comprises all data in the regime covered by BaBar, 



i.e., [1 40] GeV^. As a result, we obtained the confi- 
dential regions of different characteristics of the pion DA 
(Gegenbauer coefficients, derivatives of ipTr{x) at a; = 0, 
and its moments) by fitting the experimental data within 
the framework of LCSRs. The predictions obtained from 
the CELLO, CLEO, and the BaBar data up to 9 GcV^ 
are in good agreement with the previous fits, based only 
the CLEO data [H, HI, giving preference to an 

cndpoint-supprcssed pion DA [g. Beyond 9 GeV^, the 
best fit requires a sizeable coefficient og that inevitably 
leads to an endpoint-enhanced pion DA. The data anal- 
ysis tells us that the inclusion of the high-Q^ tail of the 
BaBar data affects mainly the Gegenbauer coefficient ag, 
while 02 and change only insignificantly. The good 
description of the experimental data up to 9 GcV^ using 
LCSRs becomes considerably less accurate after the in- 
clusion of the high-Q^ data but yields an acceptable value 
of Xndf ~ 1- This effect has been discussed before in [s^l 
at a qualitative level. The results obtained with the in- 
clusion of the high-Q'^ tail of the BaBar data indicate a 
possible discrepancy between the result of the BaBar ex- 
periment and the method of LCSRs. Indeed, the high-Q^ 
BaBar data require a pion DA with a sizeable number of 
higher Gegenbauer coefficients a^, or alternative theoret- 
ical schemes outside the standard QCD factorization ap- 
proach, see, e.g., [3l|-[3^. Similar conclusions were also 
drawn in [37| using Dyson-Schwinger equations and in 
the recent works [3^, based on AdS/QCD. 



VI. ACKNOWLEDGMENTS 

We would like to thank Simon Eidelman, Andrei 
Kataev, and Dmitri Naumov for stimulating discussions 
and useful remarks. A. P. wishes to thank the Min- 
istry of Education and Science of the Russian Federation 
( "Development of Scientific Potential in Higher Schools" 
projects No. 2.2.1.1/12360 and No. 2.1.1/10683). This 
work was supported in part by the Heisenberg-Landau 
Program under Grant 2011, the Russian Foundation for 
Fundamental Research (Grant No. 11-01-00182), and the 
BRFBR-JINR Cooperation Program under contract No. 
FlOD-002. 



[1] S.J. Brodsky, G.P. Lepage, Adv. Ser. Direct. High Energy 

Phys. 5 (1989) 93. 
[2] A.V. Efremov, A.V. Radyushkin, Phys. Lett. B 94 (1980) 

245. G.P. Lepage, S.J. Brodsky, Phys. Rev. D 22 (1980) 

2157. 

[3] F. del Aguila, M.K. Chase, Nucl. Phys. B 193 (1981) 517; 

E. Braaten, Phys. Rev. D 28 (1983) 524. E.P. Kadant- 

seva, S.V. Mikhailov, A.V. Radyushkin, Sov. J. Nucl. 

Phys. 44 (1986) 326. 
[4] B. Melic, D. Miiller, K. Passek-Kumericki, Phys. Rev. D 

68 (2003) 014013. 
[5] A.V. Radyushkin, Dubna preprint P2-10717, 1977 



hep-ph/0410276 



[6] V.L. Chernyak, A.R. Zhitnitsky, Phys. Rept. 112 (1984) 
173. 

[7] A.P. Bakulev, S.V. Mikhailov, N.G. Stefanis, Phys. Lett. 

B 508 (2001) 279; Erratum; ibtd. B590, 309 (2006). 
[8] H.J. Behrend et al, Z. Phys. C 49 (1991) 401. 
[9] J. Gronberg et al, Phys. Rev. D 57 (1998) 33. 
[10] B. Aubert et al, Phys. Rev. D 80 (2009) 052002. 
[11] A. Khodjamirian, Eur. Phys. J. C 6 (1999) 477. 
[12] K. Nakamura et al, J. Phys. G 37 (2010) 075021. 
[13] A.P. Bakulev, S.V. Mikhailov, A.V. Pimikov, N.G. Ste- 
fanis, Phys. Rev. D 84 (2011) 034014. 



N.G. 



Stefanis, Phys 
Stefanis, Phys. 
Stefanis, Phys 
Phys. C 44 (1989) 157, 



Rev. 



N.G. 



N.G. 



Lett. 



Rev. 



[14] S.V. Mikhailov, N.G. Stefanis, Nucl. Phys. B 821 (2009) 

291; corrected in arXiv:0905.4004v6. 
[15] S.S. Agaev, V.M. Braun, N. Offen, F.A. Porkert, Phys. 

Rev. D 83 (2011) 054020. 
[16] A.P. Bakulev, S.V. Mikhailov, 

D 67 (2003) 074012. 
[17] A.P. Bakulev, S.V. Mikhailov, 

B 578 (2004) 91. 
[18] A.P. Bakulev, S.V. Mikhailov, 

D 73 (2006) 056002. 
[19] V.M. Braun, I.E. Filyanov, Z. 
[20] S.S. Agaev, Phys. Rev. D 72 (2005) 114010. 
[21] P.d. A. Sanchez, it et al., 1101.1142. 
[22] V.M. Braun et al, Phys. Rev. D 74 (2006) 074501. 
[23] R. Arthur et al, Phys. Rev. D 83 (2011) 074505. 
[24] P. KroU, Eur. Phys. J. C 71 (2011) 1623. 
[25] A. Schmedding, O. Yakovlev, Phys. Rev. D 62 (2000) 

116002. 

[26] S.J. Brodsky, G.F. de Teramond, Phys. Rev. D 77 (2008) 
056007. 



[27 

[28; 

[29' 
[30 

[31 
[32' 
[33^ 
[34 

[35' 
[36 

[37 
[38; 



S.V. Mikhailov, A.V. Pimikov, N.G. Stefanis, Phys. Rev. 
D 82 (2010) 054020. 

N.G. Stefanis, Nucl. Phys. Proc. Suppl. 181-182 (2008) 
199. 

M.A. DonneUan et al, PoS LAT2007 (2007) 369. 

S. V. Mikhailov and N. G. Stefanis, Mod. Phys. Lett. 

A24, 2858 (2009). 

A.V. Radyushkin, Phys. Rev. D 80 (2009) 094009. 
M.V. Polyakov, JETP Lett. 90 (2009) 228. 
A.E. Dorokhov, Phys. Part. Nucl. Lett. 7 (2010) 229. 
Y.N. Klopot, A.G. Oganesian, O.V. Teryaev, Phys. Lett. 
B 695 (2011) 130. 

A. Stoffers, L Zahed, arXiv:1104.2081 

X.-G. Wu and T. Huang, Phys. Rev. D82, 034024 (2010); 

larXiv:1106.4365] [hep-ph]. 

H.L.L. Roberts et al Phys. Rev. C 82 (2010) 065202. 
S.J. Brodsky, F.-G. Cao, G.F. de Ter amond, Phys. Rev. 
D 84 (2011) 033001; [mXiv:1105.3999I 



